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ABSTRACT 


The  electron  beam  fluorescence  technique  was  used  for  performing 
spatially  resolved  density  measurements  of  a  ternary  gas  mixture  com¬ 
posed  of  helium,  nitrogen,  and  hydrogen.  The  flow  field  studied  was 
produced  by  a  multinozzle  array,  and  total  gas  density  values  exceeded 
1017  molecules/cc  within  certain  regions  of  the  flow.  The  experimen¬ 
tal  apparatus,  the  procedures,  and  the  method  of  data  reduction  used 
are  discussed. 
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NOMENCLATURE 


A 

Esju 


C'[Vij>] 


d 

EB 

I 

KS; u 

kxY^i) 


Angstrom  ( 1 0~  ®  cm) 

Coefficient  for  s*h  filter,  gaseous  specie  as  defined 
by  Eq.  (6a) 

Excitation  coefficient  for  s**1  filter,  ju^-h  specie,  including 
quenching  effects 

Wavelength- dependent  spectral  sensitivity  factor  as  de¬ 
fined  by  Eq.  (1) 

Nitrogen- helium  nozzle  throat  diameter 

Electron  beam  energy 

Electron  beam  current 

Excitation  parameter  defined  by  Eq.  (8) 

Quenching- rate  constant  of  i**1  state  of  specie  X  by 
specie  Y 
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kpp'  W 

L 

N 

Nj(l-) 
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nP 

V 

V' 

n^r) 

^ss' 

— ► 
rli 

?2j 

SX<ij> 


T 

T' 

W(s) 

W^spec 

x 

x(i) 

* 

X 

X 

£ 

y 

A 

y 

z 

z(i) 

* 

z 


Quenching-rate  constant  of  ith  state  of  specie  p  by  specie 
p  for  radiation  observed  in  s^*1  filter 

Beam  length 
Number  density 

First  negative  electronic  system  of  N2  molecule 
0— *-1  vibrational  band  of  N^tl-) 

Number  density  of  ith  level  of  the  Xth  specie 
Number  density  of  pth  specie 

Number  density  of  ptk  specie  obtained  from  filter  signal 

Number  density  of  p*k  specie  obtained  from  spectrometer 
signal 

Position- dependent  number  density  of  p^  specie 

Ratio  of  output  signals  s  and  s' 

ith  location  of  filter  data  acquisition 
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Photon  emission  rate  per  unit  electron  beam  length 
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Gas  temperature 
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stl:l  filter  signal  divided  by  electron  beam  current 
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ith  iteration  of  x 
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Axial  flow  stream  distance 
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Value  of  z  determined  using  filter  signals 
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4 

a 

ex(y) 

as? 

As) 

Ax(ij) 

TXU) 


Coordinate  defined  by  Fig.  1 

Indicial  subscript  denoting  specie,  as  are  ju  and  ju' 
Branching  factor  of  i  to  j  transition  of  X**1  specie 

Detector  solid  angle  subtended  at  the  source 
(s) 

Product  kv  (i)  r  (i) 

HV  ju 

Wavelength  of  i  to  j  transition  of  specie  X 

Excitation  cross  section  for  g  to  i  transition  of  specie 

X  for  electron  beam  energy  Eg;  for  convenience 

cr  (gi;E  )  is  shortened  to  a  (gi) 

A  13  A 

Lifetime  of  level  i  of  specie  X 
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SECTION  I 
INTRODUCTION 


Electron  beam  fluorescence  diagnostics  of  gas  dynamic  expansion 
flow  fields  have  been  demonstrated  during  the  previous  decade  to  be  a 
valuable  source  of  information  under  low- density  flow  conditions.  Pre¬ 
vious  applications  of  the  technique  have  included  spatially  resolved 
measurements  of  gas  density  of  single- specie  flow  fields,  nitrogen  ro¬ 
tational  and  vibrational  temperatures  in  a  variety  of  flow  fields,  and 
gas  densities  of  binary  monatomic  mixtures  expanding  from  a  sonic  or¬ 
ifice  (Refs.  1-3).  Recent  applications  of  the  technique  at  AEDC  have 
extended  the  utilization  of  the  technique  to  studies  of  the  interactions 
of  two  flow  fields  (Ref.  4),  each  composed  of  a  single  gas  specie, 
and  measurements  of  both  gas  density  and  temperature  have  been  per¬ 
formed.  However,  all  of  the  aforementioned  applications  share  the 
characteristics  of  being  low-density  and  either  single  or  binary  gas 
specie  expansions.  Extension  of  the  technique  simultaneously  to  higher 
density  (i.  e.  ,  total-density  values  on  the  order  of  10^  molecules/cc) 
and  multiple- specie  flow  fields  results  in  a  substantial  increase  in  com¬ 
plexity  of  the  measurement  and  interpretation  of  the  data.  Such  an  in¬ 
crease  in  complexity  is,  of  course,  a  result  of  the  requirement  to  con¬ 
sider  collision  quenching  effects  and  the  existence  of  radiation  "cross¬ 
talk'1  through  either  optical  filters  or  spectrometers.  The  problems 
are  further  complicated  if  the  observed  radiation  system  of  one  or  more 
of  the  species  is  weak  relative  to  the  others.  In  this  event,  quenching 
constants,  background  correction,  cross-talk,  and  both  radiative-  and 
collision- induced  cascading  effects  must  be  accurately  measured.  De¬ 
spite  the  existence  of  such  difficulties,  the  requirement  of  spatially  re¬ 
solved  density  data  for  certain  testing  applications  justifies  the  utiliza¬ 
tion  of  the  electron  beam  technique.  Included  in  these  applications  are 
diagnostics  of  exhaust  plumes,  which  are  typically  multi- specie  flow 
fields,  and  studies  of  the  gas  dynamics  of  multi- specie  interacting  flows 
such  as  are  used  for  continuous  chemical  laser  systems.  It  is  for  the 
latter  application  that  this  work  was  done. 

The  gas  dynamic  source  of  this  investigation  was  an  array  of  noz¬ 
zles  and  jets  as  shown  in  Fig.  1,  Appendix  I.  The  smaller  orifices 
were  hydrogen  (H2)  gas  sources,  and  the  larger  orifices,  or  nozzles, 
produced  supersonic  flows  of  nitrogen  (N2)/ helium  (He)  mixtures,  which 
were  used  for  simulation  of  gases  actually  employed  in  the  chemical 
laser  system.  A  detailed  description  of  the  nozzle  bank,  operating  con¬ 
ditions,  and  simulation  parameters  is  given  in  Ref.  5.  Using  the  source 
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parameters  listed  in  Ref.  5,  it  was  estimated  that  the  order  of  magni¬ 
tude  of  the  free-stream  total  number  density  would  be  10^  atoms  and 
molecules/ cc.  For  this  flow  field,  it  was  desired  to  study  the  inter¬ 
actions  of  the  gas  dynamic  expansions  of  the  H2  orifices  and  the  N2/He 
nozzles.  Specifically,  it  was  desired  to  measure  the  axial  and  radial 
spatial  profiles  of  the  absolute  value  of  the  number  density  of  the  indi¬ 
vidual  species,  H2,  He,  and  N2.  It  was  also  desired  to  obtain,  if  pos¬ 
sible,  information  concerning  the  spatial  variation  of  the  gas  tempera¬ 
ture. 


The  following  sections  of  this  report  discuss  the  theoretical  and 
experimental  details  of  the  application  of  the  electron  beam  fluores¬ 
cence  technique  to  the  nonroutine  problem  at  hand  and  the  degree  of 
success  attained  in  the  effort. 


SECTION  II 
THEORY 


2.1  SINGLE-SPECIE  EXCITATION 

The  interaction  of  a  continuous,  high-energy  electron  beam  of 
energy  Eg  and  current  I  with  atoms  or  molecules  in  ground  state  g 
produces  a  steady- state  number  density  in  excited  quantum  state  i, 
n^U),  and  a  steady- state  detected  photon  emission  rate  for  beam 
length  L  for  the  i-to-j  transition,  S^Uj).  It  is  well  known  that  S^Uj) 
can  be  written  as 


C'[Ax(ij)]/3x(ij)CTx(gi;Eg)  X  I  x  n x( g)  x  L 

1  +  kxxb)rx(i)nx(g) 


(1) 


where  /3xUj)  is  the  branching  factor  of  the  i-to-j  transition,  rxd)  Is 
the  lifetime  of  the  state  i,  and  kxx^)  the  total  quenching- rate  con¬ 
stant  for  state  i.  It  is  noted  that  both  i  and  j  are  merely  two  of  a  mani¬ 
fold  of  quantum  states  for  the  specie  X  and  its  ions.  C'[A-x(ij)]  denotes 
the  effect  of  the  wavelength- dependent  optics  transmission  factor  and 
photocathode  sensitivity.  The  factor  ££2/4tt  represents  the  source  solid 
angle  subtended  by  the  detector  and  is  the  location- dependent  factor  of 
the  equation.  If  radiative  cascading  effects  are  present,  it  is  easily 
shown  that  the  only  change  required  for  Eq.  (1)  is  the  replacement  of 
ox(gi;Eg),  henceforth  denoted  by  crx(gi)>  by 
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gi)  -  o x(gi)  +  S  «7x(gt))Sx(ti)  (2) 

where  t  represents  the  manifold  of  states  for  which  the  energy  eigen¬ 
values  exceed  the  energy  eigenvalue  of  state  i.  The  t-to-i  branching 
factor  is  denoted  by  Inclusion  of  effects  due  to  both  collisional 

quenching  of  the  cascading  radiation  and  collisional  induced  cascading 
can  be  shown  to  require  only  the  replacement  of  crx(gi)  in  Eq.  (1)  by  the 
expression 


fiy(ti)  +  k  Y  Y<ii)r ..( t)n  Y(g) 

r  v(f?>)  -  ffv(Ri)  +  £  ffv(fit)  — - — - - - - - 

XB  XB  ,  L  1  +  kxx(l)rx(t)nx(p) 


(3) 


where  kxx^^  is  partial  quenching- rate  constant  denoting  a  colli¬ 
sional  transfer  process  whereby  specie  X  in  state  t  undergoes  a  tran¬ 
sition  to  state  i.  Only  exoergic  processes  are  considered,  and  the  re¬ 
strictions  on  the  energy  eigenvalues  remain  as  defined  for  Eq.  (2). 

The  inclusion  of  cascading  effects  in  Eq.  (1)  shows  clearly  the  de¬ 
sirability  of  selecting  for  use  as  flow-field  diagnostics  transitions  which 
are  substantially  free  of  both  collisional  and  radiative  cascading.  Un¬ 
fortunately,  only  the  latter  effect,  or  its  absence,  can  be  deduced  from 
existing  data,  for  collisional  quenching  constants,  both  total  and  partial, 
are  rarely  known  for  cascading  states.  Additionally,  one  must  rely  up¬ 
on  extrapolated  values  of  the  excitation  cross  sections  to  correct  for, 
or  estimate,  effects  caused  by  radiative  cascading,  and  this  is,  in  gen¬ 
eral,  possible  only  for  allowed  excitation  processes  due  to  the  paucity 
of  data  for  high-energy  excitation  cross  sections  for  first-order- 
forbidden  transitions.  Consequently,  for  diagnostics  applications  it  is 
desirable  to  select  transitions  which  are  the  result  of  allowed  excitation 
processes  so  that  radiative  cascading  effects  must  be  caused  by  excita¬ 
tion  processes  which  are  at  least  first-order  forbidden.  Subsequent 
radiative  cascading  transitions  from  states  excited  by  such  forbidden 
processes  are,  in  general,  much  weaker  in  intensity  than  radiative 
transitions  from  allowed  excitation  processes,  and,  as  a  result,  the 
second  terms  of  both  the  numerator  and  denominator  of  the  expression 
in  brackets  in  Eq.  (3)  are  ignored.  For  this  case,  one  may  write  Eq. 

(1)  as 


Sx(ij)/I 


LC'Ux{ij)]/3x(ij)|ffx(pi)  +  1  C7X( Rt)jS x< tijjn  x(k) 


1  +  kx  x(i)rx(i)nx(p) 


(4) 
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and  it  is  seen  that  S^(ij)/I,  an  experimentally  determined  quantity,  pro¬ 
vides  a  direct  measurement  of  nx(g),  provided  k^x^  TX^  *s  known. 
The  multiplicative  factor 

L  X  (l^)  x  C'[Ax(ij)]/3x(ij)  [cx(gi)  +  2<7X(gt)/3x(ti)j 
is  determined  by  an  in- situ  calibration  of  the  test  gas. 


2.2  MULTI-SPECIE  EXCITATION 


Extension  of  the  preceding  excitation  equation  for  a  ternary  gas 
mixture  of  species  X,  Y,  and  Z  provides  the  set  of  equations 


LC'[Aa(ij)]j8G<ip 

ffa(gi)  +  £  ffa(gt)j8a(ti)" 

na(g) 

1+4iV1)r«(iVg) 

(5) 


where  the  dummy  index  /u  indicates  a  sum  over  the  species.  It  is  as¬ 
sumed  that  heavy  particle  excitation  processes  are  negligible  and,  as 
before,  that  collisional  quenching  is  the  result  of  a  binary  collision. 
Additionally,  collisional  cascading  effects  are  ignored.  Equation  (5) 
shows  the  additional  complexity  due  to  mixture- quenching  effects;  now 
three  quenching- rate  constants  are  required  for  each  radiative  specie. 

For  many  applications,  as  was  the  case  for  this  test,  the  use  of 
optical  filters  for  wavelength  isolation  is  necessitated  by  requirements 
placed  upon  the  time  available  for  data  acquisition.  Typical  filters  in 
the  visible  regime  possess  bandwidths  on  the  order  of  10  to  30  A,  and 
the  use  of  such  filters  requires  that  consideration  be  given  to  optical 
cross-talk  effects.  For  the  ternary  mixture  X,  Y,  and  Z,  Fig.  2  shows 
the  likely  source  of  such  cross-talk,  namely,  electron  beam  excitation 
of  radiative  systems  of  all  three  species  with  respective  radiation  wave¬ 
lengths  within  the  filter  bandpass.  Specifically,  if  specie  a  =  X  is  to  be 
observed  using  a  specified  filter  for  the  i-to-j  transition  of  X,  Fig.  2 
shows  transitions  of  species  Y  and  Z  which  can  contribute  to  the  filter 
signal  if  and  ^(mn)  are  approximately  equal.  As  a 

result,  the  output  signal  of  the  filter  selected  for  observation  of  the 
i-to-j  transition  of  specie  a  is  the  sum  of  the  signals  resulting  from 
the  three  transitions  shown  in  Fig.  2.  Additionally,  for  the  model 
shown  in  Fig.  2  one  must  include,  in  general,  quenching  effects  on  the 
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cross-talk  signals  since  each  "impurity  cross-talk"  transition  may  suf¬ 
fer  collisional  quenching. 

Figure  3  shows  an  interaction  diagram  for  the  quenching  processes 
associated  with  the  radiating  species  of  Fig.  2.  The  parenthetical  quan¬ 
tities  denote  energy  levels  subject  only  to  the  restriction  that  the  initial 
excited  level  and  the  final  level  are  unequal.  The  specie  A(g)  represents 
any  of  the  species  X,  Y,  or  Z  in  the  ground  electronic  state.  The 
quenching  constant  for  the  process  indicated  in  Fig.  3,  k(s)(i),  is  the 

quenching  rate  of  excited  specie  X  in  level  i  for  collision  with  specie  A 
in  the  ground  electronic  state  g.  Superscript  s  denotes  that  the  i-to-j 
radiative  transition  of  specie  X  is  observed  by  filter  s. 

Assuming  three  filters  are  to  be  used  for  density  determinations 
of  the  three  species,  H2,  He,  and  N2,  one  can  denote  the  three  current- 
normalized  filter  output  signals  by  W(s),  s  =  1,  2,  and  3,  and  designate 
the  gases  H2,  He,  and  N2  by  the  subscripts  tu  =  1,  2,  and  3,  respecr 
tively.  One  can  now  write  the  filter  signals  as 


W(s) 


3 

2  k  ^s^.(i)r  (i)n 
n'=i  w  w  '  V 


]i 


»  s  =  1,  2,  3 


(6a) 


where  n^(g)  is  denoted  by  n^  and  k^^]'(i)  denotes  the  quenching- rate 
constant  of  the  i^“  state  of  the  specie  by  specie  ju'  as  observed  in 
filter  s.  The  coefficient  CSfU  is  seen,  by  comparison  with  Eq.  (5),  to 
be 


=  1 “  x  c'[VijWij)[V6i)  +  }  Vs°^(ti)]  (6b) 

and  the  value  of  (ij)]  is  dependent  on  s.  Therefore,  if  no  colli¬ 

sional  quenching  processes  exist,  then 

W(s)  =  2  C  n  ,  s  =  1,2,3  (7) 

[1=1  r 

Thus,  for  this  simplified  case  nine  constants  are  required,  jcs^j  , 
s,  n  =  1,  2,  and  3.  Also,  for  the  more  general  case  including  quench¬ 
ing  effects,  27  additional  constants  jk^j  ,  s,  h'  =  1,  2,  and  3  are 

required.  Fortunately,  in  practice  many  can  be  ignored,  or  order-of- 
magnitude  estimates  can  be  utilized. 
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As  will  be  shown  in  a  following  section,  the  set  of  equations  repre¬ 
sented  by  Eq.  (6)  can  be  solved  by  iteration,  thereby  obtaining  the  posi¬ 
tion-dependent  values  of  number  density,  n^r),  i  =  1,  2,  and  3.  Con¬ 
sideration  of  convergence  of  the  iteration  process  and  accuracy  and 
precision  of  the  data  will  be  discussed  in  the  aforementioned  section; 

It  should  be  noted  that  use  of  additional  filters  to  observe  two  or 
more  radiative  systems  of  a  given  specie  can,  in  principle,  provide  a 
desirable  feature  of  redundancy  for  the  data.  If  multiple  filters  are 
used  in  this  manner,  it  is  in  principle  possible  to  employ  modified 
least- squares  fitting  to  the  data  for  determination  of  the  set  of  number 
densities  from  the  over- specifying  set  of  filter  output  signals. 


2.3  RADIATIVE  TRANSITIONS 

As  a  result  of  the  short  lead-time  available  before  application  of 
the  technique  to  the  test,  which  is  described  in  Ref.  6,  it  was  not  pos¬ 
sible  to  perform  laboratory  measurements  to  determine  the  optimum 
radiative  transitions  and  quenching- rate  constants  required  for  the 
ternary  mixture.  Order-of- magnitude  estimates  of  unknown  quenching- 
rate  constants  and  excitation  cross  sections  were  used,  as  well  as  pre¬ 
vious  experience  with  N2  and  He,  to  select  the  radiative  transitions 
listed  in  Table  I  (Appendix  II)  for  the  density  measurements.  It  is 
noted  inJTable  I  that  two  transitions  for  He  were  selected:  4*D— ►2-'-P 
at  4922  A  and  3^P— ►2^S  at  5016  A;  two  transitions  for  H2  density 
measurements  were  selected:  Ho  at  4861  A  and  H2  G— ►B  at  4634  A; 
one  transition  was  selected  for  Nj;  density  measurements:  first 

negative  (0,  1)  band,  for  which,  along  with  the  He  50 16- A  line,  an  opti¬ 
cal  filter  was  in  hand.  Figure  4  shows  a  partial  energy  level  diagram 
for  the  three  species,  and  the  three  radiative  transitions.  He  3^P — ► 

2  Is,  N2  (1-),  and  Hjg,  which  were  used  are  indicated  thereon.  A  dis¬ 
cussion  of  the  quenching- rate  constants  and  excitation  cross  sections 
is  presented  in  Section  IV. 
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SECTION  III 

EXPERIMENTAL  APPARATUS 


3.1  GENERAL  DESCRIPTION  OF  ELECTRON  BEAM  SYSTEM 

The  system  required  for  injecting  the  high-energy  electron  beam 
into  the  test- gas  flow  consisted  of  an  electron  source,  a  power  supply, 
and  a  differential  pumping  system.  A  schematic  drawing  of  the  electron 
beam  system  installed  on  the  Thermo- Chemical  Test  Rig  (TCTR)  is 
shown  in  Fig.  5.  The  source  of  electrons  was  a  GE  Type  33  television- 
type  electron  gun,  with  an  oxide-coated  cathode,  modified  for  50-kv 
maximum  operation.  The  accelerator  grid  had  an  enlarged  opening  of 
approximately  0.  48-cm  diameter,  allowing  maximum  beam  current  op¬ 
eration.  The  power  supply  for  the  filament  consisted  of  batteries 
yielding  from  6  to  18  v  dc.  To  obtain  the  required  energy  of  40  kv,  the 
anode  of  the  gun  was  grounded,  and  the  cathode  was  connected  directly 
to  the  negative  output  of  a  0  to  50  kv,  0  to  5  ma  power  supply.  Both  the 
grid  and  the  preaccelerator  bias  voltages  were  supplied  by  batteries. 
Figure  6  shows  the  circuit  diagram  of  the  gun  power  supply. 

The  electron  gun  itself  was  mounted  within  a  leaktight,  insulated 
Lexan®housing  which  was  maintained  at  2  atmospheres  (atm)  of  SFg,  as 
was  the  canister  containing  the  filament,  preaccelerator,  and  grid  volt¬ 
age  batteries,  the  voltage  divider  network,  and  the  tygon  tube  leading 
the  power  supply  cables  from  the  battery  canister  to  the  electron  gun 
housing.  This  arrangement  prevented  arcing  and  corona  discharge. 

The  beam  was  magnetically  focused  and  deflected  to  coincide  with 
a  two- stage  orifice  opening  into  the  test  section.  The  first- stage  open¬ 
ing  into  the  test  section  was  of  0.  076-cm  diameter  and  1.  26-cm  length. 
The  second  stage  was  of  0.  127-cm  diameter  and  0.  704-cm  length  and 
was  separated  from  the  first  stage  by  0.  318  cm.  A  schematic  drawing 
of  the  two-stage  orifice  is  shown  in  Fig.  7,  and  a  photograph  of  the  or¬ 
ifice  attached  to  the  electron  beam  drift  tube  and  mounted  in  the  TCTR 
is  shown  in  Fig.  8.  As  can  be  seen  in  Fig.  8,  the  end  of  the  drift  tube 
was  water-cooled. 

The  background,  or  ambient,  gas  density  into  which  the  electron, 
beam  was  injected  typically  exceeded  1  x  1017/cc".i.  To  maintain  the 
pressure  within  the  electron  gun  at  approximately  10" 5  torr  or  less, 
as  required  for  stable  and  long-life  operation,  the  orifice  was  differ¬ 
entially  pumped.  By  referring  to  Figs.  5  and  7,  one  can  see  that  the 


7 


AEDC-TR-73-96 


second  stage  of  the  orifice  was  connected  by  the  drift  tube  (as  shown  in 
Fig.  7)  and  an  external  pumping  line  to  a  10-cm  diffusion  pump.  A 
solenoid-operated,  5- cm  gate  valve  was  used  to  isolate  this  pump  from 
the  drift  tube  when  necessary.  The  first  stage  of  the  orifice  was  con¬ 
nected  by  the  tubing  surrounding  the  drift  tube  (as  shown  in  Fig.  7)  and 
an  external  pumping  line  to  a  7.  1  liter/ sec  mechanical  pump. 

A  precise  measurement  of  beam  current  was  critical  for  an  accu¬ 
rate  density  measurement.  The  copper  collector  cup  used  for  this  pur¬ 
pose  and  shown  schematically  in  Fig.  9  was  5.  7  cm  long  and  3.  8  cm  in 
diameter.  The  interior  was  carbon-coated,  and  a  double-grid  structure 
was  located  at  the  cup  entrance.  The  grids  were  comprised  of  approxi¬ 
mately  0.  040-cm  mesh  made  of  about  0.  003-cm-diam  brass  wire  and 
separated  by  0.  318  cm.  The  entrance  diameter  of  both  grids  was  0.127 
cm.  The  interior,  or  first  grid,  was  negatively  biased  at  25  to  30  volts 
(see  Fig.  9)  to  prevent  secondary  electrons  ejected  from  the  cup  inner 
wall  from  leaving  the  cup.  The  outer  grid  and  cup  were  electrically 
connected  and  could  be  positively  biased  (see  Fig.  9).  The  entire  col¬ 
lector  cup  was  water-cooled  and  mounted  on  a  sliding  vacuum  feed¬ 
through  in  order  to  permit  adjustment  of  the  vertical  position  of  the  cup 
(see  Fig.  8). 


3.2  OPTICS  AND  ELECTRONICS 

The  optical  radiation  produced  as  a  result  of  electron  beam  excita¬ 
tion  was  gathered  as  shown  in  Fig.  10  by  a  double-convex  glass  lens 
25.  4  cm  in  focal  length  and  10.  2  cm  in  diameter  and  was  focused  onto 
the  slit  of  the  optical  system,  which  was  housed  inside  a  1.  6-m-long 
aluminum  box  mounted  on  a  table  adjustable  in  both  the  horizontal  and 
vertical  planes.  A  complete  schematic  of  the  optical  system  is  shown 
in  Fig.  11.  The  part  of  the  electron  beam  fluorescence  under  observa¬ 
tion  was  located  50.  8  cm  from  the  collection  lens,  and  the  1:1  image 
was  therefore  formed  on  the  entrance  slit  50.  8  cm  from  the  collection 
lens.  The  optical  signal  incident  on  the  slit  was  chopped  by  a  PAR® 
Model  125  chopping  wheel  at  333  Hz.  The  slit  was  aligned  with  the  long 
dimension  (10  mm)  perpendicular  to  the  beam  and  the  short  dimension 
(1  mm)  parallel  to  the  beam.  Light  passing  through  the  slit  was  gath¬ 
ered  by  a  7.  6-cm-diam,  30.  5-cm-focal-length  plano-convex  quartz  lens 
located  30.  5  cm  from  the  entrance  slit.  An  iris  diaphragm  located  in 
front  of  this  lens  allowed  the  f- number  of  the  housed  optics  system  to 
be  matched  with  the  collection  lens  f-number.  An  identical  quartz  lens 
separated  from  the  first  by  a  narrow  air  gap  focused  the  gathered  light 
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onto  the  photocathode  of  an  RCA- 7265  photomultiplier  (PM)  tube  located 
30.  5  cm  from  the  lens.  Immediately  in  front  of  the  photocathode  was 
located  a  filter  wheel  assembly  capable  of  holding  eight  filters.  One  of 
the  positions  was  filled  with  a  piece  of  frosted  glass  for  alignment  pur¬ 
poses,  and  the  other  seven  positions  were  filled  with  filters.  The  spec¬ 
tral  characteristics  of  these  filters  are  shown  in  Fig.  12.  Five  filters 
were  used  for  number  density  measurements,  and  two  filters  were 
spectrally  located  with  the  (0,  1)  band  of  the  N^U")  system  for  rotational 
temperature  determination.  The  filter  wheel  could  be  rotated  into  the 
desired  position  by  an  externally  controlled  stepping  motor. 

Power  to  the  PM  was  supplied  by  a  Fluke  405B  power  supply,,  and 
the  output  signal  of  the  PM  was  processed  by  a  PAR  HR- 8  lock-in  am¬ 
plifier  with  the  reference  signal  supplied  by  the  chopping  wheel  assem¬ 
bly.  The  PM  output  was  fed  into  a  data  logging  system  for  recording 
on  paper  and  "grocery"  tape,  or  was  fed  into  the  Y-axis  input  of  an  X- 
Y  plotter  for  obtaining  relative  number  density  profiles.  The  data  log¬ 
ging  system  also  recorded  the  PM  voltage,  lock-in  amplifier  scale, 
beam  current,  beam  voltage,  and  TCTR  parameters. 

In  order  to  more  accurately  measure  the  electron  beam  fluores¬ 
cence  of  hydrogen,  the  optical  system  in  the  aluminum  box  was  remov¬ 
ed  near  the  end  of  the  test  and  replaced  with  a  3/4-m  Spex®  spectro¬ 
meter.  A  dove  prism  was  placed  in  front  of  the  entrance  slit  to  again 
make  the  long  dimension  (20  mm)  of  the  slit  perpendicular  to  the  elec¬ 
tron  beam  image.  The  slits  were  opened  to  200  Mm,  and  the  wave¬ 
length  was  set  on  4861  A  (the  Hjj  line).  An  EMI-6256B  PM  tube  was 
mounted  at  the  exit  slit,  and  the  PM  was  maintained  at  267°K  by  cold 
gaseous  nitrogen.  The  PM  signal  was  processed  by  an  Ortec®  photon 
counting  system  which  included  an  amplifier,  discriminator,  and  dual 
counter-timer. 


3.3  FLOW  VISUALIZATION 

• 

For  flow  visualization  photographs,  a  camera  was  mounted  to  ob¬ 
serve  the  flow  field  through  a  7.  6-cm-diam  quartz  window  directly  a- 
cross  from  the  window  through  which  spectral  measurements  were  ob¬ 
tained.  Since  the  electron  beam  was  stationary  and  directed  downward 
from  the  end  of  the  two- stage  orifice,  the  flow  visualization  camera 
had  to  be  moved  synchronously  with  the  nozzle  bank  in  order  to  obtain 
an  axial  photograph  of  the  flow  field.  .This  was  accomplished  by  attach¬ 
ing  the  camera  support  rigidly  to  the  x  (axial)  nozzle  bank  drive 
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component,  and  therefore  the  camera  was  moved  simultaneously  with 
the  nozzle  bank.  Photographs  were  obtained  by  making  time  exposures 
while  traversing  a  distance  which  was  usually  about  5  cm.  Required 
exposures  for  obtaining  quality  photographs  were  accomplished  by  cam¬ 
era  f-stop  adjustments  and  by  varying  the  x- component  rate  of  travel. 
Black-and-white  photographs  were  taken  using  a  rate  of  travel  of  10.  6 
sec /cm,  and  color  photographs  were  taken  using  a  slower  rate  of  32.  3 
sec/cm. 

The  camera  used  was  a  Super  Speed  Graphic®  with  a  135-mm 
Graflex  Optar®  f / 4.  7  lens  and  a  Graphex-  Wollensak®  shutter  assembly. 

Sometimes  filters  were  used  with  data-point  photographs  to  pre¬ 
vent  the  predominant  radiation  of  some  molecular  species  excited  by 
the  electron  beam  from  exposing  the  film.  Corning  sharp-cut  3-70  and 
3-72  filters  and  a  Corning  7-59  filter  were  used,  and  the  transmission 
curves  are  shown  in  Fig.  12. 


SECTION  IV 

CALIBRATION  PROCEDURES 


4.1  EXCITATION  AND  QUENCHING  CROSS  SECTIONS  AND 
CROSS-TALK  CALIBRATIONS 


Referring  to  Eq.  (6b)  and  recalling  that  the  coefficients  [cS)UJ  ,  s.  At 
1,  2,  and  3  contain  a  common  factor  which  is  characteristic  of  the  beam 
length  observed  and  solid  angle  subtended  by  the  optics  at  the  excitation 
point,  it  is  seen  that  the  set  of  coefficients  may  be  written  as 


C 


AOL 

in 


x  T'K 


Sfi 


S,  (I  =  1,  2,  3 


(8) 


where  KSiU  now  contains  only  dependencies  upon  wavelength,  excitation 
cross  sections,  and  cascading  processes.  The  multiplicative  quantity 
T',  of  course,  is  the  quantity  which  exhibits  variation  due  to  changes 
in  both  the  PM  tube  gain  and  AS7/ 47r,  the  geometrical  solid  angle  sub¬ 
tended  by  the  detector;  the  latter  source  of  variation,  that  caused  by 
solid- angle  changes,  is  usually  the  more  significant  of  the  two,  and 
optical  realignment  is  routinely  required  during  data  acquisition  as  a 
result  of  inadvertent  collisions  with  the  optical  apparatus.  It  is  now 
obvious  that  for  the  purpose  of  calibration  it  is  required  only  to  obtain 
laboratory  measurements  of  the  ratios  of  |cs^J,  s  =  ju  =  1,  2,  and  3, 
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using  the  intensity  of  one  of  the  more  prominent  lines  for  normalization. 
Consequently,  the  in- situ  calibration  is  required  to  provide  only  the 
value  of  AQ,  once  again  using  the  radiative  transition  which  was  used 
for  normalization  of  the  laboratory  data. 

It  is  to  be  recalled  that  the  subscript  n  =  1,  2,  and  3  refers  to  H2, 

He,  and  N2,  respectively,  and  that  W(l),  W(2),  and  W(3)  refer  to  the 
H2  filter  (the  line  at  4861  A),  the  5016-A  line  filter  of  He,  and  the 
(0,1)  band  of  N^d-)  for  N2,  respectively.  The  signals  W(4)  and  W(5), 
defined  and  used  in  a  following  section,  are  omitted  merely  for  brevity. 
W(s)iab  an<^  VV(s)test  designate  laboratory  and  test  chamber  measure¬ 
ments,  respectively.  Also,  for  determining  the  ratios  of  the  coefficients 
CSjU,  the  501 6- A  line  He  filter  with  He  gas,  C22*  will  be  used  for  nor¬ 
malization.  Therefore,  the  set  of  ratios  jcSiU  /C22L  s,  y.  =  1,  2,  and  3 
must  be  determined.  For  this  purpose,  tne  laboratory  apparatus  de¬ 
scribed  in  Ref.  6  was  used  to  provide  pure- gas  samples  of  H2,  He,  and 
N2.  The  optical  filter  signals  were  determined  as  a  function  of  the  den¬ 
sity  of  the  individual  species,  and  Fig.  13  shows  a  typical  result  for  the 
excitation  of  He  as  observed  by  both  He  filters.  Using  Eq.  (6a)  for  a 
single  specie,  it  may  be  seen  that  this  equation  can  be  combined  with 
Eq.  (9)  to  give 


y*H.b  =  l(4ff/Afi) ial)  X  (l'T'xKJKl  +  k^O^OnJ 


(9) 


Therefore,  plotting  the  left-hand  side  of  Eq.  (9)  versus  n^  yields  not 
only  CSiU  =  (An/47r)iab  •  T'  ■  KSjU  from  the  intercept,  but  also 
T^(i)  from  the  ratio  of  the  slope  and  the  intercept.  Figure  14  shows 
a  typical  experimental  result.  Figure  15  shows  the  variation  of  filter 
signals  W(l)  and  W(3)  with  He  density  and,  of  course,  demonstrates 
the  existence  of  H2  and  N2  filter  signals  when  observing  only  pure  He, 
an  optical  cross-talk  effect.  Figure  16  shows  the  variation  of  all  three 
designated  filter  signals  with  H2  density,  and  Fig.  17  shows  the  varia¬ 
tion  of  the  He  filter  signal  with  N2  density.  Using  such  data  and  least- 
squares  fitting  to  the  functional  form  of  Eq.  (9),  the  set  of  ratios  JcsAi  / 
C22j>  s,  n  =  1,  2,  and  3  was  determined  and  is  listed  in  Table  II.  It 
should  be  noted  that  the  experimental  uncertainty  of  the  coefficients 
CSjU  was  less  than  approximately  ±5  percent  in  all  cases. 

The  pure-gas  quenching  rates  kxx(i)  Tx(i)  were  determined  at 

room  temperature  in  the  manner  indicated  previously.  A  complete  de¬ 
scription  of  the  experimental  apparatus,  precautions,  and  data  reduc¬ 
tion  is  given  by  Lewis  and  Williams  (Ref.  6)  for  the  study  of  the 
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quenching  of  atomic  hydrogen  by  H2.  Table  III  lists  the  quenching  rates 
employed  for  this  study,  and  those  quenching  rates  which  have  not  been 
experimentally  determined  but  have  been  estimated  based  on  either  pre¬ 
liminary  data  or,  hopefully,  educated  guesses,  have  been  denoted. 


4.2  IN-SITU  CALIBRATION 

Using  the  electron  beam  apparatus  of  the  TCTR  as  described  in  a 
previous  section,  measurements  were  made  of  the  variation  of  the  W(2) 
signal  with  He  density,  and  Fig.  18  shows  the  result.  Using  the  same 
analysis  procedure  as  described  in  the  previous  section,  the  coefficient 
C22  (test)  was  determined.  Therefore,  using  the  equality 

0^(1030/022(^80  =  C  (lab)/C22(Iab) 

the  set  of  coefficients  |cSJu{test)|,  s,  m  =  1,  2,  and  3  was  generated. 
Table  II  lists  the  results  for  this  set  of  coefficients.  Experimental  un¬ 
certainty  of  this  set  of  coefficients  is  of  the  same  order  of  magnitude 
as  that  of  the  set  of  coefficients  determined  in  the  laboratory  measure¬ 
ments.  It  should  be  noted  that  linearity  of  the  photon  emission  rate 
SxUj)  for  each  specie  with  electron  beam  current  was  verified  over  the 
range  of  densities  observed  in  the  laboratory  calibration. 

It  was  observed  during  the  in- situ  calibration  that  the  nozzle  bank 
position  could  effectively  change  the  solid  angle  of  the  observation 
optics  (Afi).  Investigation  showed  that  the  ^-position  of  the  cooled  noz¬ 
zle  bank  had  little  effect  on  the  observation  solid  angle.  However,  at 
small  values  of  x,  the  z-position  definitely  affected  the  observation 
solid  angle.  With  the  test  section  of  the  TCTR  containing  a  few  hun¬ 
dred  microns  Hg  of  pure  helium,  the  W(2)  signal  was  measured  with 
the  nozzle  bank  clear  of  the  field  of  view.  As  the  nozzle  bank  was 
moved  in  the  x- direction  step-wise  closer  to  the  electron  beam, 
axis  scans  were  made,  and  the  W(2)  signal  was  recorded.  The  3)- 
position  was  set  at  the  zero  reference.  In  this  manner,  correction 
factors  to  Afi  without  nozzle  bank  interference  were  determined  for 
the  x-positions  of  the  nozzle  bank  at  various  locations. 
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SECTION  V 

DATA  ACQUISITION  AND  ANALYSIS 


5.1  DATA  ACQUISITION  PROCEDURES 

The  location  of  the  pitot  probe  with  respect  to  the  nozzle  bank  was 
determined  at  the  beginning  of  each  day  of  data-taking.  The  axial  loca¬ 
tion  (&- location)  was  determined  by  using  a  transit  theodolite  located 
outside  the  TCTR  (see  Fig.  10)  and  focused  on  the  pitot  probe  through 
the  flow-visualization  side  window.  By  arranging  the  line  of  sight  from 
the  theodolite  to  the  end  of  the  pitot  probe  perpendicular  to  the  x-axis, 
one  could  move  the  nozzle  bank  close  to  the  pitot  probe  and  record  the 
^-position  of  the  scanner  in  order  to  provide  a  reference  x- location. 

The  ^-position  could  be  located  in  this  manner  within  about  ±0.  1  mm. 

The  vertical  and  horizontal  (y-  and  £-)  locations  were  determined  by 
making  y-  and  ^-pitot-pressure  profiles  near  the  exit  of  the  nozzle  bank, 
without  B.2,  and  reference  and  ^-locations  were  determined  by  flow 
symmetry.  The  accuracy  in  the  reference  y-  and  ^-locations  was  better 
than  ±0.  08  mm. 

The  electron  beam  was  then  injected  into  the  flow,  and  the  ^-location 
of  the  beam  was  determined  by  the  ^-position  of  the  scanner  recorded 
when  the  top  edge  of  the  nozzle  bank  moved  into  the  beam  and  decreased 
the  beam  current  to  half  value.  The  nozzle  bank  was  then  moved  back, 
and  the  horizontal  cross  hair  of  the  theodolite  was  adjusted  to  lie  along 
the  axis  of  the  pitot  probe,  whereas  the  vertical  cross  hair  lay  along 
the  electron  beam  fluorescence.  A  slight  adjustment  of  the  theodolite 
focus  brought  the  entrance  slit  of  the  optics  housing  (or  spectrometer) 
on  the  opposite  side  of  the  TCTR  into  focus.  The  vertical  and  horizon¬ 
tal  positions  of  the  optics  housing  were  then  adjusted  until  the  center  of 
the  slit  was  at  the  intersection  of  the  cross  hairs,  with  the  long  dimen¬ 
sion  of  the  slit  parallel  to  the  pitot  probe  (or  horizontal  cross  hair). 

By  this  method,  the  ^-position  of  the  observed  portion  of  the  election 
beam  was  the  same  as  the  pitot  probe  ^-position.  The  horizontal  (£) 
location  of  the  observed  portion  of  the  electron  beam  was  determined 
by  making  relative  number  density  profiles  (PM  output  recorded  on  an 
X-Y  plotter)  near  the  exit  of  the  nozzle  bank,  without  H2,  and  by  locat¬ 
ing  the  z-position  by  symmetry. 

After  optical  alignment  was  achieved,  flow  conditions  were  set,  and 
number  density  data  were  taken  at  the  desired  positions  in  the  flow  field 
using  the  data  logging  system.  At  each  position,  the  output  of  the  PM 


13 


AEDC-TR-7  3-96 


was  recorded  with  each  of  the  seven  filters  in  the  filter  wheel  in  front 
of  the  photocathode. 

For  flow  visualization  data  acquisition,  the  nozzle  bank  was  moved 
into  the  electron  beam  in  such  a  way  that  the  beam  struck  the  top  of  the 
nozzle  bank.  The  camera  shutter  was  opened  just  before  the  traverse 
began,  thereby  photographing  the  flow  as  close  as  possible  to  the  noz¬ 
zle  exit. 


5.2  ITERATION  PROCEDURE 

The  data  reduction  process  utilized  filters  1,  3,  and  4  and  the 
spectrometer  H2  data,  designated  as  signal  5.  Filter  2,  the  5016-A 
He  filter,  was  not  used,  for  reasons  to  be  explained  in  the  following 
section. 


For  the  unknown  specie  distributions  of  the  ternary  mixture,  the 
appropriate  equation  for  determination  of  the  density  is,  of  course,  Eq. 
(6a),  which  can  be  written  as 


W(s)  =  2  C'  n  ,  s  =  1,  3, 4 

fi=i  SP  #* 


where 


k(sl(i)  * 

Mi  1  A 


V(i)ry 


] 


do) 


(ii) 


Defining  the  quantities  x  and  z  to  be  nj/n2  and  n3/n4,  respectively, 
one  can  write  the  set  of  equations  (10)  as 

W(4)  =  C'42n2[l  +  (C'41/C'42)x  +  (C'43/C'42)2]  (12) 

R41  =  W(4)/W(l)  =  [C'41x  +  C  42  +  C  ^gzl/lC'jjx  +  C'12  +  C'13z]  (13) 

and 

R34  =  W(3)/W(4)  =  [C'31x  +  C  32  +  C'33z]/[C'41x  +  C  42  +  C  43z]  (14) 
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The  primary  advantage  to  be  gained  by  ratioing  filter  signals  and 
obtaining  density  ratios  x  and  z  is  the  minimizing  of  the  effects  of 
quenching  of  hydrogen  and  helium  cross-talk  on  the  H2  density  deter¬ 
mination.  After  investigation  of  several  approaches  for  solution  of  the 
three  nonlinear  equations  in  three  unknowns,  the  following  procedure 
was  adopted. 

Using  the  excitation  and  quenching  coefficients  listed  in  Tables  II 
and  III,  respectively,  a  computer  calculation  was  performed  of  filter 
output  signals  for  the  following  ranges  of  number  density: 

1013  <  n2  <  3  x  1017  atoms/cc 

0.0  <  x  <  20.0 
0.0  <  7.  <  10.0 


The  results  of  the  calculation  provide  information  regarding  the 
sensitivity  of  the  measurement  over  the  density  domain  of  interest,  in¬ 
dicate  the  effects  of  cross-talk,  and  provide  aid  in  selecting  the  method 
of  solution  of  the  set  of  equations.  Figure  19  shows  the  calculated  var¬ 
iation  of  W(4)  with  both  helium  density  and  x,  the  ratio  of  hydrogen- to- 
helium  density.  The  variation  of  W(4)  with  z,  the  ratio  of  nitrogen- to- 
helium  density,  was  found  to  be  negligible,  and  z  =  0.  10  was  chosen  for 
the  calculation  of  Fig.  19.  Figure  20  shows  the  variation  of  R41  with 
x,  with  n2  as  a  parameter.  The  ratio  R41  was  found  to  be  relatively  in¬ 
sensitive  to  an  order- of- magnitude  variation  in  z.  However,  the  sensi¬ 
tivity  of  R41  to  He  density  is  severe,  and,  as  a  result,  relatively  accu¬ 
rate  values'of  n2  must  be  used  in  the  final  determination  of  x  from  R41. 
Figure  21  shows  the  relative  constancy  of  R41  versus  n2  for  He  density 
less  than  approximately  10*6  atoms/cc  and  for  x  <0.  10.  In  addition, 
the  ratio  R34  was  found  to  be  relatively  insensitive  to  the  magnitude  of 
x  for  helium  density  values  less  than  5  x  10*6  cc"l,  as  Fig.  22  shows. 

The  method  of  solution  of  Eqs.  (12)  through  (14)  was  selected  to  be 
successive  iterations,  and  the  results  mentioned  in  the  preceding  para¬ 
graph  indicated  the  general  solution  to  follow  the  procedure  of  determi¬ 
nation  of  n2«  then  z,  and  finally,  x.  However,  because  of  cross-talk 
effects,  it  is  necessary  to  determine  whether  the  N2  and  He  signals, 
W(3)  and  W(4),  respectively,  are  caused  primarily  by  a  large  H2  con¬ 
centration  at  the  point  in  the  flow  field  being  measured.  If  all  filter 
signals  are  caused  predominantly  by  H2,  then  the  procedure  of  solution 
should  begin  with  W(l)  rather  than  with  an  attempt  to  measure  n^  by 
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using  a  cross-talk  signal  W(4).  The  actual  value  of  x  at  which  H2  be¬ 
comes  predominant  was  determined  by  consideration  of  the  equation 

W(l)  =  C'un1[l  +  (C'12/C'11)  x  (1/x)  +  (C'13/C'H)  x  (z/x)]  (15) 

Equation  (15)  was  used  to  show  that,  for  2  =  0.  10,  if  x  ^15,  the  H2  sig¬ 
nal  becomes  the  most  important,  which  defines  the  value  of  x  at  which 
the  procedure  of  solution  was  altered.  If  x  as  provided  by  R41  is  found 
to  be  greater  than  15,  the  procedure  of  solution  requires  determination 
of  nj_,  followed  by  determination  of  x,  and  finally,  z. 

Referring  to  Eqs.  (11)  and  (12),  one  recalls  that  the  density- 
dependence  of  the  results  is  attributable  not  only  to  the  explicit  depend¬ 
ence  of  Eq.  (12)  on  n2,  x,  and  z,  but  also  to  the  effects  of  quenching, 
as  Eq.  (11)  shows.  The  iterative  method  by  which  the  quenching  cor¬ 
rections  to  the  coefficients  jCg^j  was  achieved  is  shown  in  Fig.  23  by 
means  of  the  flow  chart  detailing  the  computer- calculated  iterative 
scheme.  From  Fig.  23,  it  may  be  seen  that  the  first  iteration  pro¬ 
ceeds  by  obtaining  an  estimate  of  x  from  R^,  assuming  z  =  0.  10,  and 
using  the  low-density  He  curve  of  Fig.  20.  With  this  value  of  x,  z  is 
then  calculated,  and  these  values  of  x  and  z  are  used  for  a  determina¬ 
tion  of  n2  if  x  <  15.  At  this  point,  quenching  corrections  are  applied 
to  all  cross-talk  terms  and  z  is  again  determined,  yielding  a  value  of 
113.  The  cross-talk  terms  are  again  corrected,  and  x  is  calculated 
using  R41,  giving  nj.  The  cross-talk  terms,  as  well  as  C42,  CJi,  and 
C33,  are  corrected  for  quenching,  and  the  following  iterations  proceed 
in  like  manner.  If  x  >  15,  Eq.  (15)  is  used  to  give  nj;  R41,  to  give  nzi 
and  R43,  to  give  n3.  The  cross-talk  terms  are  adjusted  for  quenching 
as  just  described,  following  each  density  determination.  The  coefficients 
C^i,  C33,  and  C42  are  adjusted  following  determination  of  the  triad 
(nij  n2,  and  n3).  Therefore,  the  procedure  may  be  summarized  by 
stating  that  the  cross-talk  coefficients  are  maintained  current  in  den¬ 
sity  estimates  as  the  iteration  proceeds,  and  the  main  filter  coefficients, 
CJl*  C42,  and  C33  are  one  step  behind  in  the  sense  that  the  i^  iteration 
for  density  employs  the  (i-l)th  values  of  density  for  quenching  adjust¬ 
ments  in  Cjj,  C42,  and  C33. 

Employing  this  iteration  scheme,  the  iterated  density  values  are 
compared  with  actual  density  values  for  which  the  "theoretical"  calcu¬ 
lations  were  performed.  Figure  24  shows  the  dependence  on  n2  of  the 
second  iteration  of  z,  z(2).  It  is  seen  that  z(2)  is  never  more  than  ten 
percent  larger  than  z  for  n2  ^5  x  10*6  atoms/cc  and  0.  01  <  x  <  1.  0. 
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Consequently,  the  iteration  procedure  is  considered  satisfactory  for  z, 
not  only  because  it  differs  by  ten  percent  or  less  over  the  0  to  1  torr 
region  at  room  temperature,  but  also  because  its  difference  from  the 
theoretical  value  is  regular  and,  therefore,  capable  of  graphical  cor¬ 
rection.  Figure  25  shows  the  variation  of  x(2),  the  second  iteration  of 
x,  with  x,  assuming  z  =  0.  10,  for  various  values  of  n2-  It  is  seen  that 
at  the  higher  values  of  He  density,  x(2)  deviates  significantly  from  x. 
Figure  26  shows  the  second  iteration  of  He  density  as  a  function  of  n2 
for  x  =  0.  10  and  10.  0,  assuming  z  =  0.  10. 

As  described  in  a  previous  section,  detection  of  the  H^  line  of  H2 
was  also  performed  using  a  spectrometer,  rather  than  a  broad  bandpass 
filter,  and  photon- counting  electronics.  The  spectrometer,  output  sig¬ 
nal,  after  being  normalized  by  the  electron  beam  current,  can  be  writ¬ 
ten  as 


W(5)  =  (16) 

where  is  defined  by  Eq.  (11).  The  coefficients  C51  and  C52  obtained 
by  an  in- situ  calibration  are  listed  in  Table  II.  It  is  seen  that  the  cross¬ 
talk  effect  has  been  reduced  by  approximately  a  factor  of  four  by  using 
spectrometric  rather  than  filter  detection.  Figure  27  shows  the  varia¬ 
tion  of  W(5)  with  x  for  various  values  of  He  density. 

One  now  has  filter  signals  for  He  and  N2  densities,  a  spectrometer 
signal  W(5)  for  H2  density,  and  a  filter  signal  W(l)  for  H2  density  which 
is  of  lesser  reliability  than  W(5).  Two  iteration  procedures  were  em¬ 
ployed  to  obtain  values  of  H2  density  using  primarily  W(5)  and  values 
of  He  and  N2  density  using  primarily  W(4)  and  W(3),  respectively. 

The  first  procedure  of  solution  is  as  follows:  From  Eqs.  (11)  and 
(16),  it  may  be  seen  that  the  set  of  density  values  (ni,  n2,  and  n3)  is 
required  for  both  cross-talk  and  quenching  corrections.  The  values  of 
nj,  n2,  and  n3  provided  by  the  filter  outputs  and  denoted  by  (n'i,  nji>, 
and  n'3)  were  used  in  the  cross-talk  term  C52  nj},  and  nj>  and  n'3  were 
used  in  the  quenching  term  of  the  denominator  of  C51.  The  value  of 
H2  density  nj  then  obtained  from  W(5)  can  be  written  as 

=  g(l  +  M/(p5]  -  g  x  (17) 
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where 


g  =  W(5)  -  c  52  x  n  2 

and 

l  _  „(5)  '  ,  „(5)  <■ 
n  -  k ]2n  2  +  K13  n  3 


Results  of  the  computational  procedure  indicated  that  exhibited 
excessively  large  sensitivity  to  the  iterated  value  of  He  density  n'2; 
therefore,  this  procedure  was  discarded  in  favor  of  the  procedure 
which  is  described  in  the  following  paragraphs. 

The  ratio  R54  of  signals  W(5)  and  W(4)  was  formed  and  is  defined 
by 


R 


54 


3 

I  C' 
11=1 


(18) 


and  Fig.  28  shows  the  variation  of  R54  with  H2  ratio  parameter,  x,  for 
various  values  of  He  density.  It  may  be  seen  from  Fig.  28  that  the 
sensitivity  of  the  determination  of  x  using  R54  decreases  with  increase 
in  He  density,  an  effect  which  is  attributable  to  both  cross-talk  and 
quenching  effects.  Defining  x"as  the  H2  ratio  parameter  obtained  from 
Eq.  (18)  and  Fig.  20,  it  is  seen  that 


*'=  [R54(C'42  +  C'43z)  -  C'52  -C'53zl/[C'51-R54C'41] 


(19) 


Using  the  computer- calculated  values  of  R41,  R34,  and  W(4),  the 
values  of  number  density  (nj,  n2,  and  n^)  were  obtained  by  the  pre¬ 
viously  described  iteration  process.  These  values  of  number  density 
were  then  employed  for  the  quenching  corrections  to  Eq.  (19)  and,  using 
the  theoretical  value  of  R54,  x"was  determined.  Figure  29  shows  the 
variation  of  x"  with  the  actual  value  of  x  for  various  values  of  n2-  It  is 
seen  that  the  discrepancy  between  x"  and  x  is  a  monotonic  function  of  x 
so  long  as  n2  ^  1  x  10^  atoms/cc.  For  larger  values  of  n2,  each. data 
point  must  be  considered  on  an  individual  basis  for  determination  of  x. 
Figure  30  shows  the  functional  behavior  of  x"(2)  with  respect  to  n2  for 
various  values  of  x,  and  the  correction-free  regime  of  x(2)  versus  n2 
is  indicated. 
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5.3  ANALYSIS  OF  TEST  DATA 

For  the  flow  field  investigated,  the  description  and  results  of  which 
are  presented  in  Ref.  5,  the  spatial  profiles  of  the  posit  ion- dependent 
specie  density  n^fr)  were  determined  using  the  apparatus  and  the  ac¬ 
quisition  and  data  reduction  techniques  described  in  previous  sections. 
As  was  previously  mentioned,  both  filter  and  spectrometer  data  were 
obtained  for  the  measurement  of  H2  density,  and  a  separate  test  entry 
was  made  to  obtain  the  spectrometric  data  of  nj.  Since  the  spatial  lo¬ 
cations  at  which  the  filter  and  spectrometer  data  were  acquired  were 
not  exactly  coincident,  it  was  necessary  to  employ  a  spatial  interpola¬ 
tion  technique  to  reduce  the  spectrometer  data.  This  was  required  be¬ 
cause  both  the  spectral  cross-talk  and  the  quenching^  correct  ion  terms 
are  present  in  the  W(5)  equation.  Therefore,  if  Mi},  i  =  1,  2,  .  .  .  N, 
represents  the  N  spatial  locations  within  a  profile  of  filter  data,  as  de¬ 
fined  by  Fig.  1,  and  |?2j},  3  =  ^  2,  . .  .  N,  the  N  spatial  locations  of 
spectrometer  H2  data  within  the  same  profile,  it  is  seen  that  the  values 
nju^2j)  are  required  to  reduce  the  spectrometer  data.  However,  only 
nl(r2j),  j  =  l,  2,  .  .  .  N,  exist.  Therefore,  it  was  necessary  to  inter¬ 
polate  the  filter  number  dens  it  v  data  n^?^),  i  =  1,  2,  . .  .  N,  to  obtain 
the  values  at  the  positions  Jr2j|,  j  =  1,  2,  .  .  .  N.  For  this  purpose,  the 

results  were  selected  so  that  \?u\<  |?2j  |<  |^li+l  I .  and,  additionally,  for 
all  cases  i  |  » |  “  rn|,  thereby  minimizing  interpolation 

error.  It  was  assumed  that  the  values  n^r^)  and  n^d^-n)  were  ac¬ 
curately  connected  by  a  linear  relationship  to  allow  prediction  of  n^Cr) 
at  any  position  between  r^  and  Having  accomplished  this,  one 

has_values  of  n^Cr^)  which  are  used  in  Eq.  (19)  for  a  determination  of 
ni(?2j),  the  spectrometer  result. 

The  initial  reduction  of  the  data  used  the  5016-A  He  filter  in  place 
of  the  4922-A  signal.  All  previously  presented  equations  and  figures 
involving  filter  4,  the  4922-A  filter,  had  corresponding  equations  and 
figures  using  filter  2.  The  iteration  procedure  used  with  W(2)  rather 
than  W(4)  was  identical  to  that  just  presented.  The  n2  results,  how¬ 
ever,  showed  an  excessively  large  value  of  helium  density  with  the  ob¬ 
vious  discrepancy  worsening  in  regions  where  the  He  density  would  be 
expected  to  be  largest.  Moreover,  the  experimental  values  of  W(2)  in 
some  regions  of  the  flow  field  were  larger  than  the  theoretical  maxi¬ 
mum  value  of  n2  which  was  predicted  by  Eq.  (10)  with  s  =  2.  Figure 
31  shows  a  typical  experimental  result  of  the  acquired  filter  signals 
W(l),  W(2),  W(4),  and  W(5)  for  Profile  3  at  x/d  =1.0  for  Test  Condi¬ 
tion  6c,  as  defined  in  Ref.  5.  The  theoretical  maximum  for  W(2)  is  on 
the  order  of  40  output  units,  and  yet  W(2)  exceeds  100  output  units  at 
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y  =  -0.  5  cm.  The  W(4)  signal  follows  quite  closely  the  qualitative  be¬ 
havior  of  W(2)  and  is  less  than  its  maximum  value.  This  latter  obser¬ 
vation  eliminates  the  possibility  of  a  beam  current  avalanche  phenom¬ 
enon  as  being  the  source  of  the  discrepancy.  The  possibility  of  cali¬ 
bration  error  is  eliminated  since  the  r\2  values  derived  from  W(2)  and 
W(4)  agree  in  the  low- density  regions  of  the  flow  field.  A  third  pos¬ 
sible  cause  of  the  anomalous  behavior  is  the  existence  of  contaminant 
radiation,  or  excessive  cross-talk.  The  He  5016-A  filter  calibration 
data  for  N2  and  H2  density  variations  did  not  exhibit  anomalously  large 
cross-talk  effects,  but  it  should  be  noted  that  the  calibration  data  were 
acquired  for  gas  number  density  values  less  than  6  x  1016  atoms/cc, 
which  is  equivalent  to  2  torr  at  300°K.  Consequently,  data  acquired 
when  n2  1  x  10^  atoms/cc  may  well  show  discrepancies  not  predicted. 
Likely  sources  of  such  cross-talk  effects  are  the  atomic  nitrogen  ion 
lines  at  5045.  1,  5025.  7,  and  5016.  4  A,  all  within  the  5016-A  filter  pass- 
band.  Because  of  these  results,  all  flow-field  data  were  analyzed  using 
the  signal  W(4). 

Figure  32  shows  the  results  of  the  density  distributions  for  Profile 
3,  x/d  =  0.  13,  and  Fig.  33,  &/ d  =  1.  0,  using  the  method  of  data  reduc¬ 
tion  just  outlined.  These  results  show  the  extremely  high  value  of  gas 
density  obtained  for  this  flow-field  profile,  and  the  regions  of  H2  flow 
are  also  clearly  seen.  The  complex  axial  and  radial  density  distribu¬ 
tions  are  presented  in  Ref.  5.  An  additional  parameter  of  interest  for 
these  studies  of  a  chemical  laser  nozzle  array  is  the  ratio,  y,  of  n^C?) 
and  the  sum  of  n2Cr)  and  n3C?).  This  ratio,  y,  is  a  measure  of  the  gas¬ 
mixing  and  nozzle  interaction  effects  between  the  flow  field  composed 
solely  of  H2  and  the  binary  mixture  flow  field  of  He  and  N2.  It  may  be 
seen  that 


>'(?)  =  x(r)/'[l  +  z(r)]  (20) 

where  the  position  dependencies  of  x,  y,  and  z  are  parenthetically  de¬ 
noted  by  ?.  The  experimental  values  of  yCr)  corresponding  to  the  data 
of  Fig.  24  are  shown  in  Fig.  34. 


5.4  ERROR  ESTIMATES 

The  confidence  level  of  these  experimental  results  is  a  complicated 
function  of  the  test  parameters  and  specie  number  density  values  of  the 
flow  field.  Consequently,  only  a  qualitative  discussion  of  sources  of 
error  is  possible. 
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The  precision  of  the  number  density  data  is  due  both  to  the  repeat¬ 
ability  of  the  spectrometric  signals  and  the  sensitivity  of  the  measure¬ 
ment.  The  filter  and  spectrometer  signals  W(s),  s  =  1,  2,  3,  4,  and 
5,  were  obtained  at  least  twice  at  each  spatial  location,  and  these  data 
were  typically  reproducible  to  within  ±5-percent  deviation  from  the 
mean.  For  the  low- density  portions  of  the  flow  field,  this  implies  a 
±5-percent  imprecision  of  the  resulting  density  data.  However,  as  Fig. 

19  clearly  shows,  this  same  value  of  imprecision  yields  ever-increasing 
values  of  imprecision  in  the  density  values  as  the  specie  densities  in¬ 
crease.  Consequently,  one  must  conclude  that  the  density  data  are  im¬ 
precise  by  approximately  ±10  percent  over  the  regime  of  applicability 
of  the  technique. 

Values  of  accuracy  to  be  assigned  to  the  data  are  more  difficult  to 
determine.  It  is  felt  that  the  major  source  of  inaccuracy  of  the  result¬ 
ing  specie  density  values  is  a  result  of  both  the  imprecision  of  the 
quenching  constants,  which  is  typically  ±15  percent  for  the  most  signif¬ 
icant  rate  processes,  and  the  inaccuracy  of  the  quenching  constants  re¬ 
sulting  from  a  neglect  of  their  temperature  dependence.  Concerning  the 
temperature  dependence  of  the  rate  constants,  it  should  be  noted  that  the 
assumption  of  hard- sphere  binary  collision  processes  yields  a  yjr  de¬ 
pendence  for  the  rate  constant  and  indicates  that  as  the  expansion  pro¬ 
ceeds,  producing  lower  gas  temperatures  and  densities,  the  effects  of 
quenching  become  less  important.  However,  no  evidence  to  support  a 
hard- sphere  hypothesis  exists,  and  for  some  species,  such  as  Njj,  con¬ 
trary  evidence  exists.  Since  the  quenching  correction  at  high  density 
values  was  greater  than  100  percent,  the  overall  uncertainty  of  the  re¬ 
sults  for  these  regions  is  perhaps  no  less  than  approximately  ±25  per¬ 
cent.  For  flow-field  regions  of  less  than  approximately  10*°  cc”*,  to¬ 
tal  density  uncertainty  values  of  approximately  ±15-20  percent  may  be 
realistic.  Additionally,  it  should  be  noted  that  values  of  the  mixing 
ratio  y  may  be  of  lower  uncertainty  than  the  just- quoted  values  as  a  re¬ 
sult  of  minimizing  quenching  corrections  by  the  signal- ratio ing  procedure. 


5.5  ROTATIONAL  TEMPERATURE  MEASUREMENTS 

As  is  described  in  Section  3.  2,  two  filters,  spectrally  located  with¬ 
in  the  (0,  1)  band  of  the  Njjd-)  system,  were  used  to  obtain  data  in  an 
attempt  to  determine  the  rotational  temperature  of  N2.  Restrictions 
placed  upon  data  acquisition  time  precluded  performing  measurements 
of  the  spectrally  resolved  rotational  structure  of  the  vibrational  bands 
of  the  Nrjd-)  system.  The  filters  used  (maximum  transmission  at 


21 


AEDC-TR-73-96 


o  o  o 

4258  A  and  4278  A,  8-A  half-width)  in  principle  provided  signals  which 
when  ratioed  yielded  the  rotational  temperature  of  N2.  These  data  were 
acquired  for  the  entire  flow  field  studied.  However,  calibration  of  the 
filter  ratio  signals  required  not  only  measurements  of  N2  over  wide 
temperature  and  density  ranges  but  also  calibration  for  quenching  ef¬ 
fects  due  to  both  He  and  H2.  Unfortunately,  time  was  not  available  to 
complete  these  calibrations,  but  the  raw  experimental  data  are  avail¬ 
able  for  reduction  at  a  later  time  if  interest  warrants  such  an  effort. 


SECTION  VI 

SUMMARY  AND  CONCLUSIONS 


The  results  shown  in  Figs.  31-34  and  in  Ref.  5  demonstrate  the 
applicability  of  the  electron  beam  fluorescence  technique  to  diagnostics 
of  a  ternary  gas  mixture  of  relatively  high  densities.  The  preceding 
discussion  and  figures  clearly  indicate  the  limit  of  applicability  for  the 
flow  field  investigated.  Useful  data  were  obtained  at  total  density  values 
exceeding  1  x  1017  cc-1,  corresponding  to  3  torr  at  300°K.  Because  of 
the  small  physical  size  of  the  flow  field,  beam- spreading  was  not  severe, 
and  the  spatial  resolution  was  adequate. 

The  apparently  excessive  values  of  z  measured  in  the  flow  cores  of 
the  N2/He  flow  field  may  possibly  be  explained  for  the  high-density  por¬ 
tions  of  the  field  as  caused  by  collisional  energy  transfer  to  N2  from 
metastable  He  atoms,  resulting  in  enhanced  N2  radiation.  Although 
gas  separation  effects  can  not  be  eliminated  as  a  possible  cause,  the 
metastable  excitation  source  hypothesis  appears  to  be  the  more  prob¬ 
able  explanation.  The  use  of  additional  radiative  transitions  for 
measurement  of  N2  density  would  provide  an  indication  of  the  validity 
of  this  suggestion,  since  it  is  believed  that  resonance  excitation  trans¬ 
fer  processes  are  a  sensitive  function  of  the  amount  of  energy  being 
transferred.  Therefore,  if  indeed  it  is  the  He  2*S  state  populating  the 
E+  state  by  a  quasi- resonance  energy  transfer  process,  then 
radiative  transitions  from  different  vibrational  levels  of  the  Nj>  B^  £+ 
state  may  yield  systematically  different  values  of  N2  density.  This 
would  confirm  the  hypothesis. 

The  degree  of  optical  cross-talk  of  the  filters  and  spectrometers 
emphasizes  the  necessity  for  accurate  determinations  of  electron  excita¬ 
tion  cross  sections  of  the  various  gas  species  typically  encounted  in 
such  test  facilities.  This  is  particularly  important  for  those  species 
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for  which  the  main  radiative  transition  systems  are  weak  relative  to 
those  of  the  other  flow  constituents.  Such  was  the  case  of  H2,  the  radi¬ 
ation  from  which  was  the  result  of  a  dissociative- excitation  process. 
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APPENDIXES 

I.  ILLUSTRATIONS 

II.  TABLES 
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Fig.  2  Excitation-Emission  Diagram 
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Fig.  4  Partial  Energy  Level  Diagram  for  Helium,  Hydrogen,  and  Nitrogen 
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Fig.  5  Schematic  of  Electron  Beam  System  Installed  on  TCTR 


AEDC-TR-73-96 


AEDC-TR-73-96 


Fig.  8  Photograph  of  Drift  Tube,  Orifice,  Collector  Cup,  and  Pitot  Probe  in  TCTR 
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Fig.  9  Electron  Beam  Collector  Cup 
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Fig.  13  Helium  Signals  as  a  Function  of  Helium  Density 
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Fig.  14  Collisional  Quenching  Plot  for  He  31  P,  He  I1  S  Collisions 
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Fig.  15  Helium  Cross-Talk  Signal  of  Nitrogen  and  Hydrogen  Filters 
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Fig.  16  Variation  of  Filter  Signals  with  Hydrogen  Density 
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Fig.  17  Variation  of  Helium  Filter  Signal  with  Nitrogen  Density 


Fig.  18  In-Situ  Helium  Calibration 
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Fig.  19  Helium  Filter  Signal  as  a  Function  of  Helium  Density 
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Fig.  20  Ratio  of  Helium  and  Hydrogen  Filter  Outputs  versus  x 
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Fig.  23  Iteration  Flowchart 
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Fig.  24  Variation 
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Fig.  27  Variation  of  Spectrometer  Hydrogen  Signal  with  Hydrogen  Ratio  x 
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Fig.  29  Spectrometer  Ratio  x"  versus  x 
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TABLE  I 

RADIATIVE  SYSTEMS  AND  FILTER  SPECIFICATIONS 


^\Wavelength, 

\A 

Specie 

4278 

4634 

4861 

4922 

5016 

h2 

(1)  h2  G— MB 

(2)  25  A 

20A 

He 

41D—*,-21P 
85  A 

3XP— *^21S 
100  A 

*2 

No(0,  1) 
35A 

Legend:  (1)  Primary  transition  at  indicated  wavelength 

(2)  Filter  bandpass  evaluated  at  10-pcrcent  transmission  points 
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TABLE  II 

EXCITATION  COEFFICIENTS  AND  RATIOS 


S  \ 

1 

2 

3 

1 

(1)  1.  17  (-2) 

(2)  9.  50  (-18) 

3.  36  (-2) 

2.  72  (-17) 

1.  67  (-2) 

1.  35  (-17) 

2 

6.  54  (-3) 

5.  29  (-18) 

1.  00 

8.  09  (-16) 

6.  08  (-2) 

4.  92  (-17) 

3 

1.  23  (-3) 

9.  92  (-19) 

4.  47  (-2) 

3.  62  (-17) 

1.  99 

1.  61  (-15) 

4 

1.  15  (-2) 

9.  30  (-18) 

2.  22  (-1) 

1.  80  (-16) 

4.  38  (-2) 

3.  54  (-17) 

5 

1.  78  (-14)* 

8.  74  (-15)* 

0.  0 

U>  <WC22 

(2)  Cutest) 

and  C52  expressed  in  units  of  counts/ sec/ 
molecule/jua;  all  others,  in  counts/ sec/ 
molecule/  ma. 

Entries  are  to  be  interpreted  as  6.  54(-3)  = 

6.  54  x  IQ"3. 
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TABLE  III 

COLLISIONAL  QUENCHING  COEFFICIENTS 


s  N. 

1 

2 

3 

1 

(1)  3.  20  (-17) 

(2)  5.  0  (-18)** 

(3)  5.  0  (-18)** 

1.0  (-18) 

0.  00* 

1.0  (-18) 

1.  0  (-18) 

1.  0  (-18) 

1.  42  (-18) 

2 

6.  3  (-18)* 

1.  0  (-18)** 
1.  0  (-18)** 

5.  0  (-18)** 
1.  65  (-17) 

0.  00* 

1.  0  (-18)** 

1.  0  (-18)** 

0.  00* 

3 

5.  2  (-18)* 

1.  0  (-18)** 
1.  0  (-18)** 

1.  0  (-18)** 
0.  00* 

1.  0  (-18)** 

5.  0  (-18)** 

5.  0  (-18)** 

1.  65  (-17) 

4 

1.  0  (-17)* 

1.  0  (-18)** 
1.  0  (-18)** 

1.  0  (-18)** 
8.  0  (-18) 

1.  0  (-18)** 

1.  0  (-18)** 

1.  0  (-18)** 

1.  0  (-18)* 

5 

3.  20  (-17) 

5.  0  (-18)** 
5.  0  (-18)** 

1.  0  (-18)** 
0.  00* 

1.  0  (-18)** 

1.  0  (-18)** 

1.  0  (-18)** 
1.42  (-18)* 

Legend:  (1) 

(2) 

<3) 

e.  g..  KiV  =  3.20  (-17)  is  3.  20  x  1017 
cc  /  molecule 

*  Preliminary  experimental  measurements 
**  Estimated.  | 
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